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We report our investigation of the electronic and magnetic excitations of SrsNilrOe by resonant 
inelastic x-ray scattering at the Ir L 3 edge. The intra-t 2 g electronic transitions are analyzed using 
an atomic model, including spin-orbit coupling and trigonal distortion of the IrOe octahedron, con¬ 
fronted to ab initio quantum chemistry calculations. The Ir spin-orbital entanglement is quantified 
and its implication on the magnetic properties, in particular in inducing highly anisotropic magnetic 
interactions, is highlighted. These are included in the spin-wave model proposed to account for the 
dispersionless magnetic excitation that we observe at 90 meV. By counterbalancing the strong Ni^^ 
easy-plane anisotropy that manifests itself at high temperature, the anisotropy of the interactions 
finally leads to the remarkable easy-axis magnetism reported in this material at low temperature. 

PACS numbers: 71.70.Ej, 78.70.Ck, 75.10.Dg 


Since the discovery of spin-orbit induced Mott insula¬ 
tor in Sr 2 lr 04 [1], it has become apparent that electron 
correlation effects can be important in 5d transition metal 
oxides when combined with spin-orbit coupling (SOC), 
despite rather wide d electronic bands. In a perfect oc¬ 
tahedral environment of bd^ low-spin systems, e.g. in 
Ir^+ oxide compounds, SOC splits the t 2 g band occupied 
by the 5 electrons, into a fully-occupied jeff = | band 
and a narrow singly-occupied jeff = \ band. This band 
can be gapped by a fair amount of Hubbard repulsion 
driving the system toward an insulating state. This can 
persist even when the ideal jeff = \ state is not realized, 
i.e. when it is mixed with the jeff = | states due to the 
Ir4+ regular octahedral environment (see Fig. 1(b)), 
as it is generally the case in the known iridates [2-5]. 
The influence of SOC is not limited to transport prop¬ 
erties. Another very interesting perspective in the field 
of iridates is the strong anisotropy of the magnetic inter¬ 
actions produced by spin-orbital entanglement. This is 
expected to give rise to novel exotic magnetic states and 
excitations, still to be discovered [6, 7]. 

In this context, the compounds of the family AsMM’Oe 
(A = alkaline-earth metal, M,M’ = transition metal) 
are of interest because the M’ site can be occupied by 
Ir^+. Most members of the family crystallize in the 
K 4 CdCl 6 -derived rhombohedral structure in which al¬ 
ternating distorted trigonal prismatic MOe and nearly 
octahedral M’Oe coordinations constitute chains along 
the rhombohedral axis (see Fig. 3) arranged on a trian¬ 
gular lattice. They show complex magnetic behaviors 
attributed mainly to frustration and low-dimensionality, 
enhanced by strong magnetocrystalline anisotropy and 
the presence of two magnetic species. SrsNilrOe displays 
particularly intriguing properties: a complex antiferro¬ 


magnetic structure below a Neel temperature T^v of 70 K 
with slow spin dynamics and two consecutive correlated 
magnetic regimes when lowering the temperature [8-11]. 
The ordered magnetic moments are aligned along the c- 
axis, but there is currently no explanation for the ori¬ 
gin of the large anisotropy (single ion effect or exchange- 
driven mechanism) and its relation with the very high 
coercive field [12] of this material. The role played by 
SOC and anisotropic exchange is therefore a key issue. 
Three groups have performed ab initio electronic struc¬ 
ture calculations [13-15] confirming SrsNilrOe Mott in¬ 
sulating state. Nevertheless, they disagree about the 
nature of the frontiers bd electronic levels, the sign of 
the nearest-neighbor Ni-Ir interactions, and the magne¬ 
tocrystalline anisotropy. The role of SOC has been inves¬ 
tigated in the related Cu compound by resonant inelastic 
x-ray scattering (RIXS) [16, 17]. However, the distor¬ 
tion of the Ir^+ octahedral environment in SrsCuIrOe is 
monoclinic, which leads to nonequivalent Ir-0 distances 
and O-Ir-0 angles. The resulting jeff = | and jeff = | 
mixed states and their overlap with the Cu 3d orbitals are 
then at the origin of an unusual ferromagnetic exchange 
anisotropy arising from antiferromagnetic superexchange 
[17]. Ir^+ octahedra in SrsNilrOe, instead, only experi¬ 
ence a pure elongation along the trigonal c-axis, which 
preserves equivalent Ir-0 distances. SrsNilrOe is accord¬ 
ingly a nice candidate to study the influence of simple 
distortion on the electronic and magnetic properties. 

Hereafter, we report RIXS experiment on SrsNilrOe 
probing simultaneously its electronic and magnetic ex¬ 
citations. The local atomic effective Hamiltonian of 
SrsNilrOe is derived, complemented by quantum chem¬ 
istry ab initio calculations. The analysis of the magnetic 
excitations reveals a large Ising anisotropy in the Ni-Ir 
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FIG. 1: (Color online) (a) Incident energy dependence of the SrsNilrOe RIXS spectra at the I /3 Ir edge measured at 
300 K, showing the resonance of the t 2 g and Cg levels, (b) Splitting of the bd levels due to the octahedral 
environment (lODq), to its trigonal distortion and to SOC. (c) RIXS spectra at the reciprocal space position 
Q = (10 -1 7) between 10 and 300 K. The scans are offset by 75 counts for clarity. The dashed curves on the 300 K 
data correspond to least-square refinements of the excitations using Pearson-VII functions. The vertical lines on the 
10 K curve corresponds to the eigenvalues of H calculated for A=396 meV, A=294 meV, o:x^y=0 and az= 0.1 eV. 

(d) RIXS map measured at 10 K along the (10 -1 L) direction. The dotted lines show the spin-wave excitations 
calculated with J^x = Jyy = 20 meV, Jzz=^b meV and D=9 meV. The Ir and Ni main contributions are in white 

and grey respectively. 


magnetic exchange which competes with a sizable single¬ 
ion planar anisotropy for the Ni ions. This competition 
explains the directional crossover of the magnetization 
observed at high temperature. 

RIXS measurements were performed at the ID20 beam 
line of the European Synchrotron Radiation Facility. The 
beam line is equipped with a 2 m arm spectrometer, 
based on spherical Si(844) diced crystal analyzers. The 
combination of the Si(lll) double-crystal monochroma¬ 
tor with a channel-cut allows for a flexible choice of the 
energy-resolution down to 25 meV [18] at the Ir I /3 edge. 
RIXS has proven an ideal tool in the study of the iri- 
date electronic structure [3, 16, 19, 20], as it directly 
probes bd states via two successive dipole transitions 
[ 21 , 22 ]. Measurements were performed between 300 and 
10 K on a single crystal of Sr 3 NiIr 06 , grown using the 
flux method [ 11 ]. Fig. 1 (a) shows a RIXS map obtained 
from inelastic scans at different incident photon energies. 
The resonance of the Ir^+ electronic levels is observed 
at 11.214 keV for the t 2 g and at 11.218 keV for the Cg. 
For all measurements presented below, the incident pho¬ 
ton energy was fixed to 11.214 keV, which enhances the 
excitations within the t 2 g manifold. 

The RIXS spectra of Sr 3 NiIr 06 measured at the mo¬ 
mentum transfer Q = (10-1 7) of a forbidden Bragg 
peak position are shown Fig. 1(c). We first focus on the 
room temperature measurement where four features are 
visible beside the elastic peak, at the energies of 50(5), 
322(20), 568(2) and 728(5) meV, none of them showing 
any detectable dispersion within the lowest instrumen¬ 


tal resolution. The most intense peaks at 568(2) and 
728(5) meV are interpreted as d-d excitations within the 
t 2 g levels [4, 16, 17], i.e. as transitions of the hole from 
the ground state to the two lower lying filled states (see 
Fig. 1 (b)). In addition to the SOC splitting of the ^ 2 ^, 
a trigonal crystal field indeed further splits the jeff = | 
ground manifold. This effect can be quantified using a 
local atomic model for the single hole in the ^T 2 g states 
with the Hamiltonian: 

^ = A(^K)(aig|-i|ef)(ef|)+AL-S- aA 

iEx,y,z 

( 1 ) 

where A describes the t 2 g orbitals splitting due to the 
IrOe octahedron trigonal distortion (A > 0 for an axial 
elongation) and A is the SOC parameter. The symmetry- 
adapted \aig) and |e^^) wave-functions for trigonal sym¬ 
metry are defined in the Sup. Mat. [23]. The lODq 
splitting between t 2 g and Cg levels is of the order of 4 eV 
(see Fig. 1(a)), i.e. much larger than A and A, and there¬ 
fore the eg levels can safely be ignored. The last terms, 
aiSi = —2Jii < S'X* > Si {i = x,y, z), are the molecular 
field components produced on the Ir spin by its two Ni^+ 
nearest neighbors. For the analysis of the 300 K data, far 
above T^v, we first assume ax,y,z = 0 , i.e. weak influence 
of the magnetic interactions. 

The values of A and A are obtained by constraining the 
eigenvalues of H to the energies of the t 2 g manifold exci¬ 
tations determined from the experiment. Two solutions 
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exist, depending on the sign of A: A=396(l) meV and 
A=294(7) meV for the first one and A=417(4) meV and 
A=-218(8) meV for the second one. The wave-functions 
(shown in Fig. 2) for the highest doublet in the sub¬ 
space (occupied by the hole) are: 

In d) + |e'g 4) 

< . ( 2 ) 

|Q _ p|Qig4)-|eg+,t)-4eg 4) 

with the evolution of p with A/A shown in Fig. 2. The 
sign of A can not be determined by the RIXS analysis 
alone. However, the Ir magnetic moment refined in neu¬ 
tron diffraction [ 11 ] seems in better agreement with the 2 : 
component of the total magnetic moment, < Lz + 2Sz >, 
calculated for a positive A (0.25 ps) than for a nega¬ 
tive one (1.35 /i^), as shown in Fig. 2. Nevertheless, a 
strong reduction of the magnetic moment by quantum 
fluctuations, enhanced by low dimensionality or geomet¬ 
ric frustration, is not excluded. 

Another argument in favor of A > 0 comes from ab ini¬ 
tio calculations. For these, we used a quantum chemistry 
inherited computational scheme [3, 4, 16, 24], in which it 
was shown that the combination of Multireference Con¬ 
figuration Interaction method with SOC (MRCI+SOC) 
is able to reproduce the RIXS spectra. The local Ir-5d 
electronic structure has been investigated in the basis of 
multi-configurational wave-function-based methods [23] 
which is implemented in MOLCAS 7.8 code [25] to en¬ 
sure a proper description of the multiplet physics. 

The 54 spin-orbit Ni-Ir coupled states and their en¬ 
ergies were calculated for the cluster shown in Fig. 3 
[23]. In order to estimate the excitations only within 
the Ir-5d orbitals, we performed calculations by replac¬ 
ing the Ni^+ with non-magnetic Zn^+. Without SOC, 
we observe three doublet states, in which the degenerate 
‘^Eg doublets lie 0.21 eV higher than the ground state 
(see Table I). This is in perfect agreement with the en¬ 
ergy windows calculated with Ni^+, confirming that these 




No SOC 

With SOC 

Exp. values 

Configuration 

eo ei 

£2 

(M 

0 

eo ei 62 



0.00 0.21 

0.21 

0.00 0.64 0.80 

0.00 0.57 0.73 

IV’i) 

^ 2 ) 

V’s) 

= d\d\d^ 
— d\d^d\ 

1.00 

1.00 

1.00 

0.46 0.54 0.00 
0.27 0.23 0.50 
0.27 0.23 0.50 

0.58 0.42 0.00 
0.21 0.29 0.50 
0.21 0.29 0.50 


TABLE I: SD-MRCI results for the t 2 g level splittings 
(in eV) and weights of the different Ir-t 2 ^ 
configurations. In the absence of SOC, ei = 62 
corresponds to the energy of the doubly degenerate Eg 
state. The energies and weights extracted from the 
experiment for A >0 are shown in the last column. 


correspond to the excitation energies within the Ir- 5 t 2 ^ 
manifold. Furthermore, the ground doublet state corre¬ 
sponds to the ^{^ 2^3 configuration, where di is the aig 
orbital that can be written as dz 2 , while ^2 and d^ cor¬ 
respond to the Eg orbitals in the trigonal representation. 
One can hence conclude that the aig orbital lies higher 
in energy than the e'g ones, in accordance with a positive 
A. 

Taking into account the SOC, two excitations are cal¬ 
culated at 0.64 and 0.80 eV. These values are remarkably 
close to the energies found in the RIXS spectra, confirm¬ 
ing that the associated peaks correspond to the splitting 
of the Jeff = I quartet into two doublets. The weights 
of different t^g contributions to the spin-orbit coupled 
doublet states are shown in Table I. The ground state is 
not an equal mixture of the t 2 g orbitals, the weight of 
l^pi) being larger than the other two degenerate config¬ 
urations, with a ratio of 0.46 : 0.27 : 0.27 (to be com¬ 
pared with 0.33 : 0.33 : 0.33 in a perfect cubic environ¬ 
ment). Overall, the room temperature RIXS measure¬ 
ments confronted to the ab initio calculations are fully 
consistent and allow us to describe quantitatively the 
electronic spectrum of the Ir^+ and its departure from the 
perfect jeff = | state due to its non-cubic environment. 
An Ir- 5 t 2 ^ orbitals splitting A ranging between 0.21 and 
0.30 eV is obtained, leading to a ratio A/A « 0.53 — 0.74 
(A=0.396 eV estimated from RIXS). This scheme does 
not explain the weak signal observed at 322 meV. Its 
origin is unclear, but a feature with similar energy and 
spectral weight was reported in almost all the iridates 
[3, 4, 16], suggesting that it might be an intrinsic char¬ 
acteristics of this class of materials. 

Lowering the temperature down to 10 K, the two 
intense peaks related to the t 2 g electronic transitions 
change shape and acquire a substructure (see Fig. 1(c)). 
This reflects the increasing influence, as the tempera¬ 
ture decreases, of the Ni-Ir nearest neighbor magnetic 
exchange interactions on the Ir electronic states. The 
molecular field produces a splitting of the peaks that can 
be calculated from the last term of Eq. 1 . Although the 
features observed in the range 500-900 meV do not allow 
to univocally determine the values of cXx,y,z^ an exam¬ 
ple of splitting is shown on the 10 K spectrum of figure 
1(c), in reasonable agreement with the experiment. It is 
obtained with ax^y=0 and az=0.1 eV, agreeing with the 
strong uniaxial character of the magnetic structure, and 
with the value of the Jzz interaction determined here¬ 
after. 

At 300 K, the excitation at the lowest energy is broad 
with a position of the maximum at ~50(5) meV. It 
rises in intensity while narrowing as the temperature is 
lowered (see Eig. 1(c)), and its energy rapidly shifts to 
?^90(5) meV at 100 K, below which it remains constant. 
This excitation cannot be accounted for by the local elec¬ 
tronic level scheme derived previously. It is attributed to 
a magnetic excitation whose observation by RIXS is al- 
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FIG. 2: (Color online) Evolution of the ^ components of 
the orbital <L^>, spin <Sz>^ total magnetic moment 
< Lz ^ 2Sz > and of p as a function of y. The grey 
stripe shows the Ir magnetic moment refined from 
neutron diffraction [11]. The amplitudes of the 
probability of presence of the t 2 g hole are shown at the 
bottom: for a perfect octahedral environment (y =0), 
for a trigonal compression (y = —0.52) and for a 
trigonal elongation (y = 0.73) of the octahedron. The 
latter two are in agreement with the RIXS results but 
only the positive y is compatible with quantum 
chemistry calculations. The surfaces of constant 
amplitude of probability are colored from red for pure 
\aig^^) state, to blue for pure state, 

through white for an equal mix of these states. 


lowed in this system [23]. The absence of any detectable 
dispersion is illustrated by the RIXS measurement along 
the (10 -1 L) direction (see Fig. 1(d)). Another mag¬ 
netic excitation was identified at 35 meV by inelastic 
neutron scattering (INS) on a SrsNilrOe powder sam¬ 
ple [26]. It was also shown to persist up to 200 K, i.e 
much above Tjv, a behavior attributed to the presence of 
strong spin correlations associated to the low-dimensional 
nature of the magnetism. We succeeded in reproducing 
the spin-wave modes at 35 and 90 meV observed by INS 
and RIXS respectively using the Holstein-Primakov for¬ 
malism within the linear approximation (see Fig. 1(d)) 
[23]. We used a simple model, in line with the reported 
magnetic structure [11, 23], of isolated chains with (i) 
strong anisotropic antiferromagnetic exchange interac¬ 
tions between the nearest-neighbor Ni and Ir ions along 
the chain: J^^=20(2) meV, Jzz=^^{‘^) meV and (ii) 

magnetocrystalline anisotropy of Ni^+: DS^ term with 
D=9{1) meV corresponding to an easy-plane single-ion 
anisotropy perpendicular to the chain axis. The effective 
interchain interaction was estimated from both the order¬ 
ing temperature and the coercive field. It is more than 
one order of magnitude smaller than the intrachain in¬ 
teractions, affecting negligibly the spin excitations. It is 
furthermore irrelevant to the spin wave dispersion along 


the chains. It can thus be safely ignored [23]. Extract¬ 
ing the ion-dependent neutron spectral weight of the ex¬ 
citations clearly demonstrates that Ir(Ni) mainly con¬ 
tributes to the high(low) mode [23]. This explains the ab¬ 
sence of the INS 35 meV mode in the RIXS data since 
the measurements at the Ir edge are mainly sensi¬ 
tive to the resonant inelastic signal of Ir. Our model is 
highly constrained by the energies and dispersion of the 
modes extracted from the RIXS and INS data [23] and 
can uniquely account for the experimental results if one 
considers this easy-plane single-ion term. This term is 
a hallmark difference between the physics of SrsNilrOe 
and SrsCuIrOe, which otherwise shows rather similar 
anisotropic exchanges [17]. The strength of the exchange 
anisotropy is driven by p, which is related to the rela¬ 
tive weight of the \aig) and |e^) wave functions in Eq.2. 
Those are respectively preserving or not Lz when flipping 
Sz, hence contributing to the isotropic (resp. anisotropic) 
part of the interactions [17]. The parameter p varies with 
the octahedral distortion A and might be tunable by a 
chemical or external pressure. 



FIG. 3: (Color online) Magnetic susceptibility of 
SrsNilrOe calculated by Monte Carlo simulation using 
the spin-wave parameters. The blue (red) curve is 
obtained for the magnetic field applied along 
(perpendicular to) the chain direction. Top insert: 

Comparison with the measurements reported in Ref. 

[11]. The arrows indicate the crossing of the curves. 
Bottom insert: Detail of the SrsNilrOe structure with 
the central distorted IrOe octahedron along the chain 
axis. 

To further support this model, it is found that the 
Ni anisotropy perfectly matches the one reported for 
the isostructural compound SrsNiPtOe, for which D « 
9 meV [27] and which is even considered as a prototypal 
example of a "large-D” system. Furthermore, the com¬ 
petition between strong exchange anisotropy {Jzz/Jxx) 
and single-ion anisotropy {DS^) on the Ni site produces 
a hallmark feature in the magnetization measurement: 
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at high temperature, the susceptibility measured for a 
maguetic field applied aloug the c-directiou is lower thau 
the oue obtaiued with a field iu-plaue. Below 200 K, the 
two curves cross due to the iucreasiug iuflueuce of the 
large exchauge auisotropy wheu correlatious build-up, ul¬ 
timately leadiug to a uuiaxial autiferromaguetic order. 
The directioual crossover iu the maguetic susceptibility 
is uicely reproduced by Moute-Carlo simulatious usiug 
a oue-dimeusioual chaiu model aud the parameters ex¬ 
tracted from the spiu-wave calculatious (see Fig. 3) [23]. 
Fiually, the large coercive field reported iu SrsNilrOe is 
simply explaiued as the miuimum field respousible for 
the cohereut reversal of the spiu chaius, its value agree- 
iug with the iuterchaiu coupliug [23]. 

To summarize, we have determiued the bd Ir elec- 
trouic scheme iu the SrsNilrOe chaiu compouud through 
our RIXS measuremeuts aud their aualysis. The re- 
sultiug spiu-orbital eutaugled grouud state produces a 
stroug auisotropy of the iuteractious, which is evideuced 
iu our aualysis of the spiu-wave excitatious iucludiug the 
two gapped modes seeu by RIXS aud ueutrou scatter- 
iug. This uuiaxial auisotropy of the iuteractious com¬ 
petes with a stroug Ni^+ siugle iou easy-plaue auisotropy, 
a uuique feature of this system iu the series. Compar- 
iug our maguetic Hamiltouiau to the oue of a receut uu- 
merical study of a mixed-spiu XXZ chaiu with siugle-iou 
auisotropy, it is fouud that, iu the parameter space of ex¬ 
chauge aud siugle-iou auisotropies, SrsNilrOe lies close to 
the bouudary separatiug the observed ordered Isiug ferri- 
maguetic phase from a disordered XY phase [28]. These 
results demoustrate the poteutial of SOC iu triggeriug 
emergiug physics. Iu mixed 3d-5d trausitiou metal chaiu 
systems, it cau allow to probe uuexplored regious of their 
phase diagram aud quautum phase trausitious. 
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Supplemental Materials: Anisotropic interactions opposing magnetocrystalline 

anisotropy in SrsNilrOg. 
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DEFINITIONS OF THE ORBITALS USED IN THE SIMPLE HAMILTONIAN ACCOUNTING FOR THE 

RIXS MEASUREMENTS 
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AB INITIO CALCULATION DETAILS 

The local li-bd electronic structure has been investigated in the basis of multi-configurational Complete Active 
Space Self-Consistent Field (CASSCF) [Rl] method which is implemented in MOLCAS 7.8 code [R2] to ensure a 
proper description of the multiplet physics. To make such multi-electron approach feasible, a small part representing 
the crystal in which a finite number of atoms are treated explicitly in the calculation, called a ‘cluster’, has to be 
chosen. In our case, we include an IrOe octahedra, its neighboring NiOe trigonal prismatic coordinations and the 
adjacent Sr atoms (see insert of Fig. 2 in the article). This choice is expected to provide an accurate description of 
the charge distribution in the region of IrOe octahedra. To model the crystalline environment, we employ a standard 
embedding procedure of the cluster [R3, R4] including (i) a set of optimized point charges [R5], and (ii) the Ab Initio 
Embedding Model Potentials [R6] to avoid artificial excessive polarization of the cluster electrons towards the point 
charges. An active orbital space has been chosen to include the t2^-like orbitals of the li-bd shell and the magnetic 
orbitals of the Ni atoms to describe the excitations within the Ir-5t2^ manifold. The chosen active space is enlarged 
to include also the Ir-5e^ orbitals when the t 2 g -T Cg excitations are to be considered. The dynamical correlation is 
then taken into account at a Single and Double Multireference Configuration Interaction (SD-MRCI) level using the 
Restricted Active Space Self-Consistent Field [R7] module of MOLCAS, in which a maximum number of 2 electrons 
are allowed to be excited from the 2p orbitals of oxygen atoms surrounding the iridium, without relaxing the orbitals 
from CASSCF. Further enlarging the external active space in the SD-MRCI step, e.g. to include all doubly occupied 
Ni-Sd shells or the empty t 2 g orbitals from the Ir-6d shell to improve the electron correlation description, does 
not significantly change the results. The scalar relativistic effects are accounted for with the Douglass-Kroll-Hess 
Hamiltonian [R8, R9]. The wavefunctions and energies of the spin-orbit free multiplets obtained at the SD-MRCI 
level are then included for the spin-orbit coupling (SOC) calculation using the Spin-Orbit Restricted Active Space 
State Interaction method [RIO, RII] as implemented in MOLCAS code. The cluster ‘molecular’ orbitals are expanded 
with ANO-RCC basis sets [RI2] in all calculations, with the contraction schemes (8s7p5d2/), (6s5p4d), (453pld), 
and (655p2dl/) for the Ir, Ni, O, and Sr atoms, respectively. 

As expected, due to a large crystal field splitting lODg of about 3.7 eV (estimated from Single and Double MRCI 
+ SOC, in good agreement with the experimental value of 4 eV), all 5 electrons of Ir-5d are arranged in a low-spin 
S = \ manner to occupy the t2^-like orbitals. In the spin-orbit free framework, these three doublet states of Ir^+ and 
the S = 1 states of the two Ni^+ ions can couple to produce three sets of one sextet, two quartet and two doublet 
states, lying in energy windows of 0 — 33 meV and 0.24 — 0.26 eV. This gives an estimate to the Ir-5t2^ splitting of 
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0.2 — 0.24 eV due to the non-cubic crystal field. Allowing these multiplets to interact via SOC results in a total of 
54 spin-orbit coupled states lying in energy windows of 0 — 22 meV, 0.64 — 0.66 eV, and 0.82 — 0.84 eV. In order 
to get better estimate on the excitations only within the Ir-5d orbitals, further calculations have been performed by 
replacing the Ni^+ with non-magnetic Zn^+ ions, as presented in the main article. One result of these SD-MRCI 
calculations is that the three lowest doublets are always dominated by t^g configurations, weighting more than 98% 
of the total contributions, confirming a minor t 2 g — Cg mixing to the ground state wavefunction. 


MAGNETIC EXCITATION AND SPIN WAVES CALCULATIONS 

We checked beforehand that single magnon excitations in SrsNilrOe can effectively be detected by RIXS. The RIXS 
scattering intensity is given by the Kramers-Heisenberg formula [R13] and is proportional to the matrix elements of 
the polarization dependent dipole operator = ef as: 


^huj^huj' (x|(5d/|C,,|3p|)(3p||Ce|5di)^ (E3) 

where 5di and 5df stand for the initial and final d states and Spa for the intermediate states. Using the dipole 
matrix reported in ref [R14], we can compute the matrix elements displayed in Eq. E3 for the different possible 
intermediate Spa states. They are gathered in Table I for initial and final d states in the ground level with respect to 
d-d excitations (z.e. |0,t) and |0,^) states defined in the main text). All of them are not null thereby implying that 
single magnons with spin deviations in the ground level do contribute to the scattering. 
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TABLE I: Computed (0,t (i)|C'e'|3p|)(3p| ICgjO, | (t)) quantities for the different 3p| states. 

To model the spin dynamics, spin-wave calculations have been performed using the Spinwave software developed at 
Laboratoire Leon Brillouin for inelastic neutron scattering experiments. Based on the Holstein-Primakov approxima¬ 
tion, the code diagonalises any bilinear spin Hamiltonian and takes into account (isotropic or anisotropic) exchange 
couplings acting between neighboring spins, as well as single ions anisotropy terms. The initial spin configuration is 
obtained from the mean field solution of the following Hamiltonian: 

n = Y^J^, (5f D{Stf (E4) 

i iENi 

which describes an isolated chain with alternating spin 1 and spin ^ antiferromagnetically coupled along the 
chains [R15]. The exchange tensor acting between the first neighbors along the c axis is highly anisotropic, 
defined by Jxx and Jzz- 'H also takes into account the magnetocrystalline anisotropy of the Ni spins, modeled 
by and which corresponds to an easy plane anisotropy perpendicular to the 2: axis. The point symmetry 

between Ni and Ir site is 3, thus allowing for antisymmetric exchange interactions with a Dzyaloshinskii-Moriya 
D-vector pointing along the c-axis. The upper limit estimate for this interaction, 5 meV, shows that it is not 
the leading term in the Hamiltonian, and has thus been neglected, together with the interchain interactions (see below). 

Eigure El shows the calculations performed with the parameters: Jxx = 20 meV, Jzz = 46 meV and D = 
9 meV, along a Q-direction corresponding to the zone explored during the RIXS experiment (see left Eig. El). Our 
calculation shows the presence of two excitations around 35 and 90 meV with a small dispersion of the order of 10 meV. 

We also calculated the powder average of the spin wave spectrum in order to compare our model with the powder 
neutron inelastic experiment performed bu Wu et al. [R16] (see right Eig. El). Because of the integration over the 
momentum directions, the resulting spectra lack directional information but give a weighted density of states for each 
magnitude of Q. This results in two band-like excitations, the lower one reproducing fairly well the INS results. The 
values of the different parameters are highly constrained by the results from the RIXS and INS experiments. Indeed 
the position of the maxima of the two modes are equal to Jzz — D and to 2Jzz^ whereas the dispersion is related to 
Jzz'- The amplitude of the dispersion decreases as Jxx decreases and as Jzz increases. 
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The Spinwave software also allows one to identify the contributions arising from the different magnetic species. 
Indeed, since a spin-wave is a collective excitation, each spin within the unit cell contributes with however a different 
weight. A close look at those individual contributions shows that the lower branch is essentially due to the Ni 
contribution while the upper branch arises mainly from the Ir contribution (see Fig. F2). This explains why we did not 
observed the low-energy excitation in the RIXS experiment as the energy was fixed at the Ir-I /3 resonance edge. The 
reason why the high-energy excitation was not seen in INS is probably due to the weakness of the signal in the rather 
high Q-range investigated. Overall, the calculations give an excellent agreement with both RIXS and INS experiments. 
Note that Figure F 2 shows two interlaced branches around 25 and 90 meV. This is due to the fact that the unit cell 
contains two Ir and two Ni atoms. Only one branch actually carries significant physical intensity as shown in Figure FI. 


Using this program, we also examined the influence of the mean interchain interactions on the spin wave spectrum 
by calculating the response of three chains on a triangular lattice: two of them are antiferromagnetically coupled, 
while the third one is unfrozen in orientation. This spin configuration is one of the two that are compatible with 
the neutron diffraction experiment [R15], the other one involving a modulation of the amplitude of the magnetic 
moments, which cannot be implemented in the Spinwave software. Note that these two magnetic structures are 
equivalent up to a global phase. The inclusion of the interchain interactions modifies slightly the energy of the 
two branches (see Fig. F3) by increasing (decreasing) their energy for ferromagnetic (antiferromagnetic) interchain 
coupling without affecting the amplitude of the dispersion. 


The interchain exchange interaction might be estimated from the ordering temperature, but this is not straightfor¬ 
ward for mixed-spin quantum chains where the intrachain exchange anisotropy competes with single-ion anisotropy 
so that the spin dimensionality evolves from almost XY to almost Ising by decreasing the temperature. An upper 
bound can nevertheless be perceived through analytical and numerical estimates from different models [R17-R21]. 
In the case of classical spin-S chains without anisotropy, by treating exactly the effects of the intrachain exchange 
J|| and in mean field those of the interchain exchange the ordering temperature is analytically computed in the 
form T/v ~ J\\S{S + 1)((8/3)^_L'>^_l/^||)^^^ where 2:_l is the number of chains with which each chain interacts. With 
S = 1 / 2 , J|| = Jxx = Jyy = 20 meV and T/v « 80 K ^ 6.9 meV this leads to the exchange ratio r] = J±/J\\ ~ 0.02 
and an interchain exchange J± ~ 0.4 meV. With the larger exchange Jzz a smaller exchange ration 77 is obtained. 
Neel ordering is inhibited by interchain correlations and quantum fluctuations, which would suggest that J_\_ is larger. 
Approached numerically these effects in the case of the extreme limit of S=l /2 quantum spins lead indeed to an 
increase of the exchange ratio up to 77 ~ 0.1 and therefore of interchain exchange up to ~ 2 meV. These decrease 
to 77 ~ 0.075 and J_\_ ~ 1.5 meV for quantum spins S=3/2. Neel ordering on the other hand should be favored by the 
reduction of the spin dimensionality caused by the single ion and intra-chain exchange anisotropies that should lead 
to much lower values of the interchain exchange J_l, as suggested by the exchange ratio 77 ^ 0.009 obtained in the 
Ising limit. It follows that the J± ~ 0.4 meV estimated from the classical spin-S chains without anisotropy can be 
considered safely as providing an upper bound. 


As a matter of fact, another independent estimate of J± can be obtained from the coercive field He which 
was observed up to 55 T in the magnetization processes [R12]. We assume that this coercive field corresponds 
to the minimum field of coherent reversal of spin chains. In case of the frustrated magnetic structure with two 
frozen chains of antiparallel moments over three with the third being paramagnetic, there is a degeneracy of the 
possible configurations. Each frozen chain can be surrounded by (3 + n) chains of opposite moment orientation and 
6-(3+n) chains of same moment orientation, with n = 0, 1, 2, 3. On the other hand, each unfrozen chain is always 
surrounded by 3 chains of opposite moment orientation and 3 chains of same moment orientation. We can write 
(3 + (n) t3)Jperp ~ C He with C ^ 0.058 for He given in Tesla and Jperp in meV, which leads to values of Jperp (~ 
1.06 meV for (n) = 0, ^ 0.80 meV for (77.) = 1, « 0.64 meV for {n) = 2 and « 0.53 meV for (n) = 3) of the order of 
the upper bound extracted from the Neel temperature. Accordingly, the interchain exchange is expected to have a 
minor effect on the spin wave spectrum and confirms the quasi-one dimensional character of SrsNilrOe- 


We also calculated the effect of the intrachain Dzyaloshinskii-Moriya interactions, whose component along the c 
axis is allowed. We found that the small changes produced by this term in the energy position and the dispersion of 
the spin waves cannot be seen in the RIXS and INS experiments, up to values of 5 meV. Up to this upper bound. It 
has also no effect on the magnetic configuration. 
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MONTE CARLO CALCULATIONS 

Monte-Carlo simulations have been performed by considering a system of chains of alternating classical spin 1/2 
and spin 1 sites. Interchain interactions have been neglected for simplicity. The exchange interactions and single ion 
anisotropies have been fixed to that provided by the spin-waves analysis. The simulations have been performed using 
a conventional Metropolis algorithm, cooling the system from 2000K down to IK. The temperature was decreases 
exponentially (0.95 times the previous temperature). At each step in temperature, 10^ spin flips per spin were 
conducted for equilibration and a further 10^ steps for data taking. The magnetic susceptibility was calculated using 
the fluctuation-dissipation theorem: 


Xa = 


(m|) - {rugf 
NT 


(E5) 


where N is the number of spins in the simulations and (a=x,y,z is the direction. Simulations using supercells of 
3x3x200 up to 3x3x800 have been used. The results presented in the main article correspond to the largest supercell 
with 20000 spins. 


MAGNETIZATION MEASUREMENTS 

The easy-plane nature of the magnetocrystalline anisotropy of the compound can also be observed in the magnetic 
field dependence of the magnetization (see Fig. F4). Indeed, at 300 K, the magnetization measured in the direction 
of the c axis is smaller than in the perpendicular direction. At 150 K, the inverse is observed. 


* Corresponding author: lefrancois@ill.fr 
^ Corresponding author: a.m.t.pradipto@gmail.com 
[RI] B. O. Roos, P. R. Taylor, and P. E. M. Siegbahn, Chem. Phys. 48, 157 (1980). 

[R2] F. Aquilante, L. De Vico, N. Ferre, G. Ghigo, P.-A. Malmqvist, P. Neogrady, T. Pedersen, M. Pitonak, M. Reiher, 
B. Roos, et al., J. Comp. Chem. 31, 224 (2010). 

[R3] A. Sadoc, C. de Graaf, and R. Broer, Phys. Rev. B 75 , 165116 (2007). 

[R4] R. Maurice, A.-M. Pradipto, C. de Graaf, and R. Broer, Phys. Rev. B 86, 024411 (2012). 

[R5] B. Roos and U. Wahlgren, MADPOT and MADFIT program (1969). 

[R6] Z. Barandiaran and L. Seijo, J. Chem. Phys. 89, 5739 (1988). 

[R7] P.-A. Malmqvist, A. Rendell, and B. O. Roos, J. Phys. Chem. 94, 5477 (1990). 

[R8] M. Douglas and N. M. Kroll, Ann. Phys. 82, 89 (1974). 

[R9] B. A. Hess, Phys. Rev. A 33, 3742 (1986). 

[RIO] P.-A. Malmqvist and B. O. Roos, Chem. Phys. Lett. 155 , 189 (1989). 

[Rll] P.-A. Malmqvist, B. O. Roos, and B. Schimmelpfennig, Chem. Phys. Lett. 357, 230 (2002). 

[R12] B. O. Roos, R. Lindh, P.-A. Malmqvist, V. Veryazov, and P.-O. Widmark, J. Phys. Chem. A 109 , 6575 (2005). 

[R13] T. Aberg and B. Crasemann, in X-ray Anomalous (Resonance) Scattering: Theory and Experiment^ edited by K. Fisher, 
G. Materlik, and C. Sparks (Elsevier, Amsterdam, 1994). 

[R14] E. M. E. de Groot, P. Kuiper, and G. A. Sawatzky, Phys. Rev. B 57 , 14584 (1998). 

[R151 E. Lefrancois, L. C. Chapon, V. Simonet, P. Leiay, D. Khalyavin, S. Rayaprol, E. V. Sampathkumaran, R. Ballou, and 
D. T. Adroja, Phys. Rev. B 90, 014408 (2014). 

[R16] W. Wu, D. T. Adroja, S. Toth, S. Rayaprol, and E. V. Sampathkumaran, ArXiv e-prints (2015), 1501.05735. 

[R17] T. Oguchi, Phys. Rev. 133, A1098 (1964). 

[R18] D. J. Scalapino, Y. Imry, and P. Pincus, Phys. Rev. B 11, 2042 (1975). 

[R19] J. Villain and J. M. Loveluck, 38, 77 (1977). 

[R20] E. Boersma, W. J. M. de Jonge, and K. Kopinga, Phys. Rev. B 23, 186 (1981). 

[R21] C. Yasuda, S. Todo, K. Hukushima, E. Alet, M. Keller, M. Troyer, and H. Takayama, Phys. Rev. Lett. 94 , 217201 (2005). 
[R22] J. Singleton, J. W. Kim, C. V. Topping, A. Hansen, E.-D. Mun, S. Ghannadzadeh, P. Goddard, X. Luo, Y. S. Oh, S.-W. 
Cheong, et al., ArXiv e-prints (2014), 1408.0758. 





5 



Q = (10,-1,L) Momentum Transfer (A 


FIG. FI: (Color online) Calculated spin waves with exchange interactions Jxx = Jyy = 20 meV and Jzz = 46 meV 
and magnetocrystalline anisotropy D = 9 meV for a single crystal (left) and for a polycrystalline sample (right). 
The white zone represents the region which was not accessible in the INS experiment from Fig.4 of ref. [R26] up to 

50 meV. 



FIG. F2: (Color online) Contributions of the Ir ions (left) and Ni ions (right) in the calculated spin waves with 
exchange interactions J^x = 20 meV and Jzz = 46 meV and magnetocrystalline anisotropy D = 9 meV. The 

colorbar is in log scale. 
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Q = (10,-1,L) Q = (10-1.L) Q = (10-1,L) 

FIG. F3: (Color online) Calculated spin waves with intrachain interactions J^x = 20 meV, Jzz = 46 meV, 
magnetocrystalline anisotropy D = 9 meV, and interchain interaction J± = -1 meV (left), 0 meV (middle) and 

+1 meV (right). 



FIG. F4: (Color online) Magnetic field dependence of the magnetization measured with the magnetic field applied 
parallel (filled symbol) and perpendicular (empty symbol) to the c-axis at 300 K (in cyan) and 150 K (in red). 










